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Partonic structure of constituent quark (orvalon) in the Next-to-Leading Order is used to calculate 
pion structure function. This is a further demonstration of the finding that the constituent quark 
structure is universal, and once it is calculated, the Structure of any hadron can be predicted 
thereafter, using a convolution method, without introducing any new free parameter. The results 
are compared with the pion structure function from ZEUS Coll. Leading Neutron Production in 
e + p collisions at HERA. We found good agreement with the experiment. A resolution for the issue 
of normalization of the experimental data is suggested. In addition, the proportionality of Ff and 
, F%, which have caused confusion in the normalization of ZEUS data is discussed and resolved. 
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£ ■ I. INTRODUCTION 

a : 

ZEUS Collaboration at HERA has recently published [1] data on pion structure function, F% , using the leading 
neutron production in e + p collision. The data do suggest that there is a simple relation between proton structure 
function, F%, and the pion structure function F£ ■ This assertion is believable and points to the direction that there 
! , exists a more basic and universal structure inside all hadrons. 

In Ref.[l] the normalization of <7 77r and hence, the pion structure function is fixed by two different methods: (a) 
dominance of one pion exchange, and (b) use of the additive quark model. The two normalizations differ by a factor 
\& • of two. The additive quark model makes no statement about the leading baryon production while the extraction of 
.FJ is completely based on meson exchange dynamics. A criticism of this procedure is given in [2]. It is evident that 
the issue of normalization of the data is more uncertain than it was thought before. ZEUS Coll. now believes that 
most likely the final result will lie between the two options [3]. The issue is what to use for F{t) that parameterizes 
the shape of the pion cloud in the proton. In this paper we will offer an alternative, which does not rely on those 
assumptions and renders support to a normalization, which lies between the two options used by the ZEUS Coll. 
*-£h \ In what follows we have used the structure function of a Constituent Quark (CQ), (in Ref. [4] it is called valon and 
hence the valon model) which is universal to all hadrons and from there, with a convolution method, the structure 
O i| function of pion is obtained. The motivation for such an approach is based on the fact that such a model for soft 
production has found phenomenological success at y/s < 100 GeV and low pr- It should, however, be understood that 
1 when s is high enough to generate a significant component of hard subprocesses (y/s > 200GeU), a tot and average 
, pt will both increase and inclusive distributions will lose their scaling behavior. Nevertheless, the soft component is 
unchanged, and the model remains valid. 

Our knowledge of hadronic structure is based on the hadron spectroscopy and the Deep Inelastic Scattering (DIS) 
■ data. In the former picture quarks are massive particles and their bound states describe the static properties of 
hadrons; while the interpretation of DIS data relies upon the QCD Lagrangian, where the hadronic structure is 
intimately connected with the presence of a large number of partons. It has been shown [5] that it is possible to 
perturbatively dress a valence QCD Lagrangian field to all order and construct a constituent quark in conformity 
with the color confinement. Therefore, the assumption of constituent quark as a valence quark with its cloud of partons 
is reasonable. The Cloud is generated by QCD dynamics. A complete description and calculational procedure for 
obtaining the constituent quark structure, and hence, the hadronic structure functions are detailed in [6]. 

II. FORMALISM 

By definition, a valon is a universal building block for every hadron. In a DIS experiment at high enough Q 2 it is 
the structure of valon that is being probed, while at low enough Q 2 this structure cannot be resolved and it behaves 
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as the valence quark and the hadron is viewed as the bound state of its valons. For a U-type valon one can write its 
structure as, 

F?(z, Q 2 ) = ^z{q% +q § ) + \z{qa + qi + q% + q § ) + ... (1) 

where all the functions on the right-hand side are the probability functions for quarks having momentum fraction z 
of a U-type valon at Q 2 . These functions are calculated in Ref. [6] in the next-to- leading order and we will not go 
into the details here. Suffices to note that the functional forms of the parton distributions in a constituent quark, (or 
valon), is as follows [6]: 

zq^ k (z,Q 2 ) = az b (l-z) c (2) 



zq m (z, Q 2 ) = az\\ - zf\l + V z + Hz ' 5 ] (3) 

The parameters a, b, c, a, etc. are functions of Q 2 and are given in the appendix of Ref. [6]. The above parameterization 
of the parton distribution in a valon is for light quarks, u and d. For heavy quarks additional phenomenological 
assumptions are needed to be made. It is known that in proton the strange quark distribution is smaller than up 
quark by a factor of 2 at some regions of x and by the time x reaches down to 10 -4 ,we have xs = xu. As for the 
charm quark content of proton, the picture is less clear. Early treatments of heavy parton distributions assumed 
that for Q 2 > ttiq, the heavy quark, Q, should be considered as massless. At the opposite extreme, the heavy quark 
has never been regarded as part of the nucleon sea but produced perturbatively through photon-gluon fusion. The 
number of flavors remain fixed, regardless of Q 2 . This treatment is used in GRV 94 parton distribution [7]. Both 
schemes have their own deficiency For example, in the former scheme, the heavy quark should not be treated as 
massless for Q 2 > ttiq and in the latter scheme one cannot incorporate large logarithms at Q 2 » rag. CTEQ [8] on 
the other hand, have used an interpolation, which produces relative features of both schemes. In this treatment, the 
heavy quarks are essentially produced by photon-gluon fusion when Q 2 m TOq and considered as massless quark when 
Q 2 » nig. In our treatment, however, since we are dealing with low x- values, we will assume zs = zu, and zc — azs, 
where a is a factor taken to be equal to the ratio of the strange and charm quark masses when Q 2 < rn 2 harm : and 
a = 1 for Q 2 > m 2 harm . Although such an undertaking is not free of ambiguity, however, that docs not change our 
qualitative arguments regarding the normalization of the data. 

Equations (1-3) completely determine the structure of a valon. The structure function of any hadron can be written 
as the convolution of the structure function of a valon with the valon distribution in the hadron: 

F 2 h (x,Q 2 ) = J2 f dyG^F^i-^ 2 ) (4) 

where summation runs over the number of CQ's in a particular hadron. F^iz, Q 2 ) denotes the structure function of 
a CQ (U, D, U, D, etc.), as given in equation (1), and Gcq (y) is the probability of finding a valon carrying momentum 

h 

fraction y of the hadron. It is independent of the nature of probe and its Q 2 value. Following [4], [6], and [9], for the 
case of pion we have: 

G CQ/P , on (y) = B[jl + \ v + l f ^ vT (5) 

where G cq_ is the U- valon distribution of in 7r + as well as the D- valon distribution in n~ . Similar expression for 

pion 

G(jQ/ pionl anti- valon distribution in a pion, is obtained by interchanging /i <-> v. In the above equation B(i,j) is Eulcr 
(3 function and Gco_(y) are non-invariant distributions satisfying the following number and momentum sum rules: 

h 

I GcQ{y)dy = l W yGcQ(y)dy = l (6) 

Jo C q Jo 

For pion, the numerical values are: \i = 0.01, v = 0.06. We have, however, tried a range of values for /i and v and 
did not find much sensitivity on F% data against this variations. In [4] [10] it is estimated that fi = v = 0., which 
is very close to the values that are used in this paper. The flatness or almost flatness of the valon distribution in 
pion is attributed to the fact that the valons are more massive than the pion, so they are tightly bound. From SU (2) 
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symmetry one should expect that ji — v. In our calculation, /i is slightly different from v, indicating a small violation 
of SU(2) symmetry. This violation is very small and the data on is not sensitive enough to make a large difference. 
Significant asymmetry is observed in proton sea and its implications are discussed in Ref. [6] in the context of the 
valon model. 

The pion has two valons (or constituent quarks), for example, ir + has a U and a D, therefore, the sum in equation 
(4) has only two terms. Parton distributions in a pion, say n + , is obtained as: 

Kai.(x,Q 2 )= f ^Gu(y)u^(*,Q 2 ) (7) 

Jx y ,ri ~ v 




q^(x,Q 2 ) = [ 1 ^Gu(y) qw (-,Q 2 )+ f^G^{y) q ^{-M 2 ) (9) 

x+ Jx y " + v Jx y «+ y 

Similar relations can be written for 7r~ and ir° . In the above equations the subscripts U and D are the two valon 
types in 7r + . Equations (7, 8, 9) along with equation (4) completes the evaluation of the pion structure function. In 
figures (1) we present F£(x,Q 2 ), as calculated above, for the fixed Q 2, s corresponding to the ZEUS data of Ref.[l]. 

III. DISCUSSION AND THE DATA 

In the previous section we outlined the procedure for calculating F£(x,Q 2 ). The results are shown in figures (1) 
by the square points. From the figures, it is evident that for smaller x and lower Q 2 values, the results of the model 
calculations are closer to the additive quark model normalization of the ZEUS Coll. data (see figure 19 of Rcf.[l]). 
As we move towards the large x values the calculated structure function decreases and gets closer to the effective 
flux normalization of the data. If we are to trust in the valon model results, which provides a very good description 
for the wealth of proton structure function data as well as other hadronic processes, we can conclude that the two 
normalizations used by ZEUS may be relevant to different kinematical regions and lends support to the assessment 
made by the ZEUS collaboration that the final results will lie between the two options [3]. It is true that the valon 
model resembles the additive quark model in that the contributions from each valon are added up. But here we 
are mainly dealing with the parton content of each valon, which is derived from the perturbative QCD. The valon 
distribution in pion serves only as a phenomenological mimic of the pion wave function, we further note that in the 
above calculation none of the ZEUS data for pion structure function is used and no data fitting is performed. 
ZEUS Coll. makes the observation that there is a simple relationship between the proton structure function and the 
effective pion flux normalization of the pion structure function (see figure 18 of Ref. [1]); namely 

F* (EF \x,Q 2 )^kFZ(x,Q 2 ) (10) 

with the proportionality constant, k = 0.361. We have calculate the right-hand side of Eq. (10) in the valon model and 
compared it with the effective flux normalization, F \x , Q 2 ) , in the left hand side. The copmarison is presented 
in figures (1). As one can see from the figures, the relationship holds rather well at all Q 2 values. To avoid any 
misleadings, we emphasizs that the direct calculation of KT (x, Q 2 ) in the valon model (Square points in figures (1)) is 
different from F2^ EF \x, Q 2 ) and hence, does not support the effective flux normalization of F%. In other words, our 
finding merely states that if we scale by a factor of k = 0.37 we arrive at equation (10). Figures (1) also indicates 
that the above relationship holds rather well at lower Q 2 values. As we move to higher Q 2 this relationship at the 
lowest x-value gets blurred, but at large x and high Q 2 it continues to hold. 

Since our model produces very good fit to the proton structure function data in a wide range of both x — [10~ 5 , 1] 
and Q 2 — [0.45, 10000] Gev 2 , we have also attempted to investigate the relation: 

FZ{x,Q 2 )*\f%{\x,Q 2 ) (11) 

which is based on color-dipole BFKL-Regee expansion and corresponds to the ZEUS's additive quark model normal- 
ization. We have calculated both sides of the equation (11) in the valon model. The results are shown by the solid 
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lines and square points in Figures (1). Although We get similar results as in Fig. 19 of Ref.[l], but this does not say 
much about the pion structure function data, because the additive quark model normalization of the data is based on 
the above equation. It only restates that our model, indeed, correctly produces proton structure function. The main 
ingredient of our model is the partonic content of the constituent quark which is calculated based on QCD dynamics. 
Convolution of this structure with the constituent quark distribution in the pion appears, to some extent, to give 
support for the additive quark model normalization of the pion structure function data. In fact, we agree with Ref. 
[3] that the final results should be somewhere between the two normalizations used by ZEUS, being much closer to 
the additive quark model scheme than the pion flux scheme. 

It is worth to note that the following relationship also holds very well between ZEUS's pion flux normalization data 
and proton structure function: 

F^ ef \x,Q 2 ) = \fV{\x,Q 2 ) (12) 

This relationship is essentially the same as Eq. (10), except for the factor | in front of cc that makes a small correction 
to the factor k of Eq. (10). In Figure (2) we have presented the right-hand sides of both equations (10) and (12) 
along with the effective flux normalization data from Ref.[l] at a typical value of Q 2 = IhGeV 2 . The same feature is 
also prevalent for the other values of Q 2 . 

Both Eqs. (11) and (12) indicate that, regardless of the choice of normalization scheme, the valence structure of 
proton, as compared to pion, is shifted to the lower a; by a factor of |, So that valence x in proton corresponds to |x 
in pion. 

ZEUS has observed that the rate of neutron production in photo-production process, in comparison to that of pp 
collision, drops to half and from this observation ZEUS Coll. has concluded that cr 77r ~ |o" 7P whereas one expects 
to get a rather |, both from Regge factorization and the additive quark modcl[2]. This poses a problem on the 
understanding of the dynamics of the interaction. A tentative resolution is that the discrepancy can be resolved if 
we suppose that in the process each valon interacts independently. That is the impulse approximation. Suppose that 
5 denotes the number of valons in the target proton that suffers a collision and Si denotes the number of collisions 
that i th valon of the projectile encounters. If we define the integer a = X^ 1 ^ then we will have 5 < a < 35. 
Furthermore, let Ps(<j) represents the probability that out of 5 independent collisions that target valons encounter, 
a valonic collisions occur. If p is the probability that cither of the other valons in the projectile also interact, then 
the probability for 5 collisions will have a binomial distribution having a — 5 valonic collision by i = 2 and 3 valons 
out of a maximum 25 possible such collisions [11]. Now, for a real photon, we can assume that it may fluctuate into 
mesons with two valons. if we denote the number of possible valonic collisions in 7 — p interaction by 5 then the 
mean number of such collisions will be 25 p whereas in pp collisions it will be 35p. The observed reduction in the rate 

of neutron production in two processes and the conclusion that <7 77r ~ implies that ^- = |. That is, there are 

twice as many collisions in pp interactions than in photo-production. 

IV. CONCLUSION 

We have demonstrated that the assumption of the existence of a basic structure in hadrons is a reasonable model 
to investigate the hadronic structure. Pion structure function measurement by ZEUS Coll. provides additional tests 
and further validation of the Model. We have presented a resolution to the issue of the normalization of the data. Our 
results are based on QCD dynamics and suggest that the correct normalization of F% is closer to the additive quark 
model normalization for low x region. The observed reduction in the rate of neutron production in photo-production as 
compared to pp collision is accounted for and concluded that there are twice less valonic collisions in photo-production 
than in pp collision. Furthermore, it appears that there are simple relations among the structure functions of hadrons; 
namely they are proportional and the proportionality ratio seems to depend on the normalization scheme chosen. 
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Figure 1 (a): Pion structure function at Q 2 = 7 GeV 2 . The circles and 
the triangles are pion flux and additive quark model normalizations 
of the data [1], respectively. Squares are the calculated results 
from the valon model ( see the text for details.) The solid line 
represents 2/3 (2/3 x, Q 2 ) and the dotted line is 
0.371 * Fp (x,Q 2 ), ooth are calculated using the model. 



1.2r 

1.1 - Q 2 = 60GeV 2 




0.1 T 

Q L I I I I I '''I I I I I I '''I 

10~ 3 10~ 2 10~ 1 

Figure 1 (d): Pion structure function at Q 2 = 60 GeV 2 . The circles and 
the triangles are pion flux and additive quark normalizations of the 
data [1], respectively. Squares are the calculated results from the 
valon model ( see the text for details.) The solid line represents 
2/3 F p 2 (2/3 x, Q 2 ) and the dotted line is 0.371 * F p 2 (x,Q 2 ), both 
are calculated using the model. 




Figure 1 (e): Pion structure function at Q 2 = 1 20 GeV 2 . The circles and 
the triangles are pion flux and additive quark normalizations of the 
data [1], respectively. Squares are the calculated results from the 
valon model ( see the text for details.) The solid line represents 



2/3 F r 2 (2/3 x, Cr) and the dotted line is 0.371 
are calculated using the model. 
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Figure 1 (b): Pion structure function at Q 2 = 15 GeV 2 . The circles and 
the triangles pion flux and additive quark model normalizationsof the 
data [1], respectively. Squares are the calculated results from the 
valon model ( see the text for details.) The solid line represents 
2/3 (2/3 x, Q 2 ) and the dottedline is 0.371 * F p 2 (x,Q 2 ), both 
are calculated using the model. 




Figure 1 (c): Pion structure function at Q 2 = 30 GeV 2 . The circles and 
the triangles are pion flux and additive quark model normalizations 
of the data [1], respectively. Squares are the calculated results 
from the valon model ( see the text for details.) The solid line 
represents 2/3 (2/3 x, Q 2 ) and the dotted line is 
0.371 * Fp (x,Q 2 ), ooth are calculated using the model. 




Figure 1 (f): Pion structure function at Q 2 = 240 GeV 2 . The circles and 
the triangles are from ZEUS Coll. [1 ] for two different normalizations. 
The squares are calculated from the valon model ( see the text for 

resent 2/3* F p 2 (2/3 x, Q ) and the dotted 
Q 2 ), both are calculated using the model. 



details.) The solid line repre 
line represents 0.371 * F 2 (x 



0.6 
0.55 E- 

0.5 - 
0.45 - 

0.4 - 



CM 



CP-35^ 

LU 

^(0.25 h 

LL 

0.2 r 
0.15 - 

0.1 - 
0.05 r 



9 







-4 



Q 2 = 15GeV 2 




J I L 



1 



10 



-3 



J I 1111 



1 



10 



-2 



Figure 2. The effective flux normalization of Pion structure function data 
as compared with the scaled proton structure function of Eqs. (1 0) and (1 2), 
The data points are from Ref . [1 ] . The solid line is 0.371 * K(x,Q 2 ) 
pertinent to Eq. (10) while the dashed line is 1/3 F p 2 (2/3 x, Q ) 
pertinent to Eq. (12), both are calculated from the valon model at Q 2 = 1 5 Ge 



